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Abstract 



We investigate the next-to-leading logarithmic (NLL) order contributions of 
the quantum chromodynamics (QCD) for exclusive processes evaluated by Monte 
Carlo methods. Ambiguities of the Monte Carlo calculation based on the leading- 
logarithmic (LL) order approximations are pointed out. To remove these ambigui- 
ties, we take into account the NLL order terms. In a model presented in this paper, 
interference contributions due to the NLL order terms are included for the gener- 
ation of the transverse momenta in initial-state parton radiations. Furthermore, a 
kinematical constraint due to parton radiation, which is also a part of the NLL or- 
der contributions, is taken into account. This method guarantees a proper phase 
space boundary for hard scattering cross sections as well as parton radiations. As 
an example, cross sections for lepton pair productions mediated by a virtual photon 
in hadron-hadron collisions are calculated, using the jet-calculus scheme for flavor 
nonsinglet quarks. 
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§1. Introduction 

Monte Carlo methods with parton shower models for initial-state parton radia- 
tions are powerful tools for the evaluation of exclusive processes for hadron-hadron 
scatterings. The parton shower models have been constructed on the basis of per- 
turbative quantum chromodynamics (QCD).I'' 

In actual calculations, the leading-logarithmic (LL) order of QCD is insufficient 
to evaluate the hadron-hadron scattering processes due to large parameter depen- 
dence, particularly the choice of a factorization scale parameter. Thus, the next-to- 
leading logarithmic (NLL) order contributionPPI should be taken into account. The 
next-to-leading order (NLO) calculation (0(as) contributions for hard processes) is 
also necessary to remove theoretical ambiguities due to the factorization procedure 
as well as the choice of the factorization scale.'S Here, is the coupling constant of 
QCD. 

Conventionally, the scaling violation of the parton distributions is calculated 
by solving the renormalization group equations in moments. Then, these solutions 
are numerically inverted to yield momentum fractions of partons. Alternatively, 
we can use parton shower models in order to evaluate the scaling violation of the 
parton distributions. One such algorithm has been proposed in Refs. [3]) andl4i. In 
this model, the scaling violation of the parton distributions is generated using only 
information from the splitting functions of the parton branching vertices and input 
distributions at a given energy. It has been found that the method reproduces the 
scaling violation of the parton distributions up to their normalizations at the NLL 
order of QCD. 

In order to allow the application of parton shower models to realistic processes, 
the matching problem should be solved, that is, double counting between the hard 
scattering and the parton showers must be avoided. So far, various methods have 
been proposed in order to solve this problem. In particular, a systematic inves- 
tigation of the factorization procedure at the NLO in Monte Carlo calculation has 
been performed, as described in Ref. [5]). In a recent work, a modified subtraction 
scheme at the NLO has been proposed.® They derived hard scattering cross sections 
with the new scheme for the Drell-Yan process as well as for the deep inelastic scat- 
tering process. Another method has been proposed in Ref. [6]), in which the collinear 
terms are subtracted from matrix elements before integration over the phase space. 
However, in these works, the initial-state parton evolutions are based on the LL 
order of QCD. 

In the conventional MS scheme,'^ it has been pointed out that matching between 
the hard scattering cross section and the initial-state parton radiation is broken in 
the calculation of exclusive processes.'^ In order to implement a factorization scheme 
that is appropriate for the evaluation of exclusive processes to the accuracy of the 

*-* In this paper, the NLL order contributions correspond to the corrections that appear from 
0(as) terms of splitting functions in parton evolutions, as well as boundary conditions due to 
momentum conservations at the collinear parton radiations. Furthermore, 0(0:3 log if^) terms of 
the decay matrix elements (three-body decay functions) for branching vertices are also included. 
Here, —K^ is a virtuality of an initial-state parton momentum. 
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NLL order of QCD, a kinematical constraint due to parton radiation is taken into 
account when considering the subtraction terms, which is cahed the MS scheme. 
The momenta of partons in the scattering processes are conserved in this method. 
As an example, cross sections for the Drell-Yan lepton-pair production mediated by 
a virtual photon in hadron-hadron collisions have been calculated.''' 

However, in the collinear region after the collinear singularity is subtracted, neg- 
ative contributions remain for the cross section obtained with both the conventional 
MS scheme and the MS scheme.® Such contributions cannot be ignored at the 
NLL order accuracy. Event generation employed in the Monte Carlo methods with a 
negative probability may not be appropriate, since the strong cancellation between 
the negative contributions from the hard scattering cross section and the positive 
contributions from parton showers may give unstable results. 

To solve this problem, the jet calculus (JC) scheme has been implemented in 
Refs. |8]) and nop . The JC algorithm has been proposed by Konishi et al. in order to 
evaluate multiparton productions due to perturbative QCD.^nJ This algorithm had 
been applied to jet evolutions as well as initial-state parton radiations. Furthermore, 
for the parton evolution in which the virtualites of partons are strongly ordered, 
additional phase space factors, compared with the MS scheme, are subtracted from 
the hard scattering cross sections.!^ 

In this paper, the NLL order contributions in initial-state parton radiations as 
well as the NLO corrections in the hard scattering cross section for exclusive processes 
using a Monte Carlo method are investigated in a consistent manner. As an example, 
we calculate Drell-Yan lepton-pair production in proton-antiproton scattering with 
the JC' scheme, where a kinematical constraint due to parton radiation for the hard 
scattering process is also taken into account. In this algorithm, most of the collinear 
contributions to the hard scattering cross section at the NLO are subtracted.® In 
order to simplify its explanation, we concentrate on the quark (g)-antiquark (q) 
annihilation process, but we include following ingredients for the evaluation of the 
physical cross sections: 

(1) Parton showers for initial-state radiations are generated at the NLL order accu- 
racy including the 0{a'^) terms for the splitting function, where the 0{al log K"^) 
contributions are also taken into account for each branching step. The scaling vio- 
lation is reproduced by evolution of the initial-state parton radiations at the NLL 
order accuracy up to the normalization of the distribution function, as well as the 
transverse momentum of partons .'^''i^ 

(2) Explicit expressions of the three-body decay functions {0{a1log K'^) terms of 
the decay matrix elements in the parton branching vertices) are included, where 
interference contributions of two gluon radiations are also taken into account in 
the exact form.f^ 

(3) Kinematical constraints due to the momentum conservation of parton radiations 
are taken into account for both factorization of mass singularity for the hard 
scattering cross section and the parton generations.'^'® 
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Since the above contributions are related to each other at the NLL order accu- 
racy, one should treat these terms consistently. Otherwise, the factorization scheme 
invariance for the physical cross sections may be broken. 

In §2, we explain our algorithm. Some numerical results are shown in §3. Section 
4 contains a summary and some comments. The explicit expressions of the formula 
implemented in our algorithm are presented in Appendices. 



§2. Basic formula for Monte Carlo algorithm 



Here, we consider the Drell-Yan lepton-pair production in quark ((;)~antiquark 
(q) annihilation in a proton-antiproton scattering 



QiPq) + QiPq) i*iq) + giPg) ^ I (p-) + i^ip+) + giPg), 



(2-1) 



mediated by a photon 7* with the virtuality (p_ +p+)^ = = Q^, where a gluon 
(g) is radiated in the final state. Here, pi {i = g, g, g) and p± denote the momenta 
of the corresponding particlesl*^ 

2.1. Subtraction scheme 

In order to obtain a finite cross section for the process qq — )• 7*5, we subtract 
collinear contributions due to branching processes q ^ qg and q ^ qg from the hard 
scattering cross section. Although these contributions are compensated by initial- 
state parton radiations, the remnant of the subtracted cross section in the NLO 
cannot be ignored at the NLL order accuracy. Here, we consider valence quarks 
(flavor nonsinglet sector) in order to avoid the quark-gluon mixing for si mp licity. 

The subtraction term divided by the Born cross section aniQ^. e r*i for the 
branching process 

<l{Pq) Q{r)+g{pg), 
with a subtraction scheme F is defined by 



(2-2) 



"'^NS 



a. 



1 



in 4 



dzd{-r'^) 27rr(l 
2e dimensions, where 



47r/i^ 



Cf 



1 - z 



(2-3) 



(24) 



and Q 



depends on the subtraction scheme 



Here, Cp = 4/3 is the color 
factor. The strong coupling constant is defined by as/i^*^ for the dimensionless cou- 



*' The Mandelstam variables are defined by s = {pq i = (p, —Pg)^, and u — {pq — Pg)"^ 

wliicli satisfy s + t + u = for massless partons. 

**' Tlie Born cross section f7o(Q^, e) is defined in Ref. [7|. 
We define 



for a function f{z) unregulated at 2 = 1. 



f{z) - 5{l - z) / dyf{y) 
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pling Us and a mass parameter ji. Here, Q|^^^(^;) = corresponds to the MS 
subtraction scheme. The momentum of the quark, r, is described in terms of 
the momentum fraction z, the virtuahty r^, and the transverse momentum as 
r = Pg — Pg = zpq + {r^/s)pg + vt with Pq ■ rx = Pq ■ vt = 0. Here, we set 
Pq — Pq — Pg — 0) because the relations — p^, —p^^p^g ^ — are expected in parton 
shower generation. 

The regularization of the infrared singularity is defined so as to conserve the 
particle number of the initial-state quarks for the branching process presented in 
Eq. (2-2). 

The parton evolution depends on the subtraction scheme for the mass singularity, 
which appears owing to collinear parton production. Including the NLL order term, 
the infrared regulated splitting function with the subtraction scheme F, denoted by 
Pi[fJ(a„z))_^, is defined bj^ 



with /3o = 11 - 2/37V/ for Nj active flavors. Here, (^i[J^^'^^^(z)j is the NLL order 

splitting function calculated with the MS scheme.^^ Rl 
The splitting function ( Pj[j^' (a^, z) j satisfies 



dz(ifs(«.,^)) =0 (2-6) 
^ 

in each branching process. Using 

f dz(P^l{a,,z)) = - f ' 
Jl-5 ^ ^+ Jo 

the nonbranching probablitiy for the nonsinglet quarks is defined by 

rK'idK^ rl-S 



dz(Pj|(a„z)) =-/ dzp[g(a„z), (2-7) 



r^2 dK' r-" , MFi, 



(2- 



Here, the partons inside an initial-state hadron have spacelike virtualities (i.e., 
kf = —Kf < 0), and 6 denotes a resolution of the momentum fraction of the final- 



state partonsl**^! The actual steps in the Monte Carlo method are similar to those 
presented in Refs. _3) andH]). 

Using this algorithm, the particle number distribution of the nonsinglet quarks 
is reproduced up to their normalization without the introduction of any nontrivial 



*•* We neglect the process q ^ q + X, since the contribution is small at the NLL order. 
**•* In Monte Carlo calculation, a cutoff parameter for small z integration is chosen, so that a 
very small s for a given hard process is not generated by parton showers. 
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weight factor. Therefore, the model guarantees the particle number for the initial- 
state valence quarks, namely. 



/' 

Jo 



dxf!l'J/^{x,K^) = N, (2-9) 



is satisfied for any K^, where x is the momentum fraction of the valence quark 
inside a proton, with Nu = 2 for the valence u quarks and = 1 for the valence d 
quark. Here, f^/pi^^ -^^) is the particle number distribution function of the valence 
quarks inside a proton calculated with a factorization scheme F. The above method 
guarantees universal (process-independent) parton distributions with a factorization 
scheme F, since the same scheme has been implemented in a deep-inelastic process.'^ 

2.2. Kinematical constraint on the hard process 

In Monte Carlo calculation, kinematical constraints must be taken into account, 
for the momentum conservation of partons in each generated event. Here, the kine- 
matical constraint due to the squared momentum of the photon = (r + pq)'^ 
to 0{as) accuracy of QCD is imposed in the delta function, in which the quantity 
f = Q'^/s is given by f = z + r^/s for on-shell gluon radiation (p^ = 0) (denoted by 
the F' scheme) as 



"•-"ns _ / "'-'ns 



dfd{—r'^) Jq dzd{—r'^ 
We define the integrated contribution as 



dz——^^5{z-f-{-r')/s). (2-10) 



d^-r"^^ 







dfd(-r2) 27rr(l-e) 



F[fs''^'^(e,f), (2-11) 



where / denotes the region of phase space for the hard scattering process being 
considered. Here, sw'^^^^^^ is a limit of the — integration. 

According to Eq. (2-10), the phase space < — < with z < 1 for the 
initial-state parton radiation covers the region < — t < (1 — 7^)5 for 1 — tm < 't < 
1(7 = S) and < — t < for f < 1 — tm{I = C) in the hard scattering cross 
section, with tm = M'^/s. The collinear contribution is also subtracted from the 
antiquark leg. Here, we consider the two regions 1 — tm < t < 1 and r < 1 — tm 
separately for the subtraction terms. 

For 1 — Tm < f < 1, we subtract the collinear contribution from the hard 
scattering cross section in the range < — t < s(l — f). If we ignore the term r^/s 
in the delta function in Eq. (2-10) and subtract the collinear contribution in the 
range < — < M^, with f = z, the matching between the hard scattering cross 
section and the initial-state radiation is broken. In this case, is larger than the 
kinematical boundary for —t{= — r^), given by s(l — f) in the region 1 — tm <1. 

2.3. Kinematical constraint on the branching processes 

The algorithm for the generation of transverse momenta in the initial-state ra- 
diation is similar to that for the Monte Carlo model presented in the LL order of 
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QCD, except that the effects due to the three-body decay functions, which are the 
O(a^) contributions, are taken into account for each branching step. 

In the flavor nonsinglet sector, the contribution of the two-gluon radiation, such 

as 



becomes large in the soft gluon region. Here, the momenta of these partons are 
denoted by k^, I1J2 and k, with if = I2 = and /cg, k'^ < 0. 

In our algorithm, the transverse momenta of partons were generated according 
to the branching probability restricted by the boundary condition determined by the 
NLL order terms. In the case of Eq. (2-12), the transverse momentum of the quark 
with momentum k is given by 



Here, we define P = (/i + 12)^ and z denotes the longitudinal momentum fraction of 
the parton momentum k for the momentum kQ. 

The upper limit of the virtuality P for the radiated partons is not properly 
determined at the LL order, since the kinematical constraint of the branching process 
is a part of the NLL order contribution. As shown in Ref|3]), the boundary condition 
is determined by the parton shower algorithm as well as the NLL order terms. 

For strong ordered virtuality of the initial-state partons, namely, —{kQ—li)'^, —{kQ — 
^2)^ < —k'^, the virtuality P is restricted by the condition P < /][Jg(— fc^) at the NLL 
order, instead of P < {1 — z)/ z{—k'^), by the kinematic boundary of the two-body 
branching, due to the fact that k^ > for -k"^ <^ —k'^,P at the LL order approx- 

imation. Here, the factor of /^g is given by the NLL terms calculated with the 
subtraction scheme F. 

In the case considered above, the two-body branching process is allowed in the 
branching steps of the initial-state parton radiations, as is usually the case in LL 
order parton shower models, except that the contributions due to the NLL order 
terms for the initial-state radiations are included in the kinematical boundary for 
the radiated partons. 

2.4. Factorized cross sections 

The hard scattering cross section with the factorization scheme F' is defined by 



q{ko) ^ gih) + gih) + q{k) 



(2-12) 



k^ = z{l - z) /Cq + 



1 - z 



(2-13) 



z 



da 



df 



-±ao{Q^O)K 



•[F'](/) 
qqiNS) 



T 



(2-14) 



with 



<7o(Q',0) 




(2-15) 



3NcQ 



Here, Nc = 3 and the electric coupling constant of the quark is defined by e^a. 
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2.4.1. Soft gluon region 

For the soft gluon region 1 — 2fM < f < 1, we consider two separate regions, 
one for 1 — fj\/ < f < 1 (/ = 51) and the other for 1 — 2tm < f < 1 — tm {I = S2). 
For I = SI, we have 

(^) = K,,{f,e) - 2F^^'P{e,f), (2-16) 

where Kqq{f, e) denotes the DY cross section integrated over < — t < (1 — t)s^^^ 
For I = S2, we have 

(^) = K^^(r,e) - 24^](^)(6,f), (2-17) 

since < s(l — f) in this region. 

2.4.2. Hard colhnear region 

The remnant of the cohinear gluon contribution of the hard scattering cross 
section integrated over the range < — t < ( or < —u < M^) for f < 
1 — 2fM (I = C) is given by 

= - (2-18) 

where -f^g^^ (t, e) denotes the DY cross section integrated over < — t < ( or 
< -n < M2).0 



§3. Numerical results in the JC scheme 



In this section, some numerical results obtained using the algorithm explained 
in the previous section are presented. In the parton shower model explained in 
the previous section, the virtuality of the initial-state partons is generated with 
strongly ordered configuration, namely, Mq < < X| < • • • , < M^. It has been 
discussed that the strongly ordered configuration of the parton virtuality specifies a 
factorization scheme,'^ which we refer to as the JC scheme. Therefore, according 
to the Monte Carlo algorithm for the initial-state radiations, we implement the JC 
scheme explained in Refs. [8]) and llOp . where we define 

QKs'(^) = l°g(l - ^) - ^^.(^) (3-1) 

in Eq. (2-3). Here, the function Pqg{z) = —Cf{1 — z) is an 0(e) term of the splitting 
function in 4 — 2e dimensions!*^! 

Since K^q (f) ^> K^-^^^-^ {f) for f < 1 — 2tm, as presented in Ref. [8]), the 
remnant of the hard collinear region after subtracting the collinear contribution 

*-* The sphtting function in 4 — 2e dimensions is defined bjJ^ 
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with the JC scheme can be safely neglected even at the NLL order accuracy. Here, 
K^g^'ir) is the DY cross section for hard gluon radiation presented in Eq. (2-1) 

integrated over < —i < {I — f)s — M^, and Kj^^j^^^^\f) is the same as that of 
the singlet sector given in Ref. 8). 

The Monte Carlo method for the Drell-Yan process that we used in this paper 
is presented in Appendix A. The explicit expressions of the factorized cross sections 
in the JC' scheme implemented in numerical calculations are shown in Appendix B. 

Here, we present the transverse momentum distributions of the virtual photon 
7* produced by a u-quark and an anti-u-quark annihilation in pp collisions with the 
center of mass energy at ^/s = 1800 GeV. 

3.1. Ambiguities at the LL order approximation 

First, we investigate the ambiguities at the LL order approximation of QCD. 

In Monte Carlo models based on the LL order of QCD, the initial-state parton 
showers are generated using the LL order splitting function, which has no compen- 
sation term for a subtraction method in the hard scattering cross section as well as 
the kinematical boundary for the hard scattering process. 

In the Monte Carlo calculation, the parton showers start from Mq = 5 GeV^ in 
both sides of protons. After the evolutions are terminated, the momenta of the gener- 
ated spacelike partons pq and pg are reconstructed. Then, the s-channel momentum 
is calculated as s = {pq +Pq)'^ = Q^- 

At the LL order of QCD, the upper limit of the virtuality for the generated 
parton is restricted by a large energy scale M, which can be chosen arbitrarily. In 
Fig. 1, the transverse momentum distributions for the virtual photon integrated over 
the range of the photon virtuality 40 GeV < Q < 60 GeV are represented by the 
solid histogram for M = 40 GeV and by the dashed histogram for M = 60 GeV with 
= Iq = 0.5 GeV^, where the two-body decay kinematics shown in Eq. (2-13) for 
initial-state parton radiations are imposed. The transverse momentum distribution 
clearly depends on the choice of M. Furthermore, the contributions only from the 
initial-state radiations are insufficient to reproduce the distribution in the hard qx 
region. Therefore, the NLO contribution should be added. 

Another ambiguity comes from the kinematical boundary in the initial-state 
parton branching vertices explained in 2.3. In Fig. 1, the dashed histogram shows 
a case with < l"^ < {I - z)/z(-:k'^) due to /cf, > in Eq. (2-13) at the LL order 
approximation for M = 40 GeVO 

In Fig. 1, the transverse momenta according to the cross section at the NLO 
contributions are also presented for M = 40 GeV by the plus symbols and for 
M = 60 GeV by the crossed symbols, which also depend on the factorization scale 
M. 

The transverse momentum distributions combined with the initial-state radia- 
tions and the NLO contributions are presented in Fig. 2. Although the M depen- 
dence of the transverse momentum distributions due to the initial-state radiations 



*■* The virtuality for outgoing partons is generated by using the nonbranching probability for 
timelike parton evolutions. 
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Fig. 1. Transverse momentum distributions for the virtual photon integrated over 40 GeV < Q < 
60 GeV with -^i = 1800 GeV. The solid and dashed histograms represent the results for 
M = 40 GeV and M = 60 GeV with P = = 0.5 GeV^, respectively. The dashed histogram 
shows a case with < < (1 - z)/z{-k'^) for M = 40 GeV. The NLO contributions are 
represented for AI = 40 GeV by the plus symbols and for M = 60 GeV by the crossed symbols. 

is partially compensated by that from the NLO contribution, the calculated results 
at the LL order approximation still depend on the factorization scale M as well as 
the phase space boundary for the branching vertices. 

3.2. The NLL order contributions 

The transverse momentum distributions including the NLL order contributions, 
explained in §2, are presented in Fig. 3. Here, the NLO contribution with the JC' 
scheme is implemented, where the phase space boundary of the hard scattering cross 
sections is imposed in subtraction terms as well as the generation of the initial-state 
parton radiations. 

Furthermore, the initial-state radiations are generated within the phase space, 

which satisfies Sm,n - Sm+l,n > {Sm,n " Sm,n+1 > -p|(„+i)) with Sm,n = 

{Pq{m) +Pq{n)f' for each branching step. Here, Pq{m){Pq{m)) ("i = 0, 1, • • • ) denotes the 
momentum of the initial-state quark (antiquark) generated in the m-th branching 
step in the forward evolution algorithm!*^ 

In Fig. 3, the transverse momentum distributions for the virtual photon in- 
tegrated over the range of the photon virtuality 40 GeV < Q < 60 GeV are pre- 

*-* The parton evolutions start from the distribution function at A/q, which is approximately 
given b}^ 

f^:^ii.,Ml) - /i^l(-,Mo^) + ^ £ f W)^ CI (f '^^o) 

with the GRV(98) distribution function for f^^^j^{x, Mi)W^ 
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Fig. 2. Transverse momentum distributions for the virtual photon integrated over 40 GeV < Q < 
60 GeV with -^5 = 1800 GeV. The sohd and dashed histograms represent the resuhs for 
M = 40 GeV and M = 60 GeV with = Iq ^ 0.5 GeV^, respectively. The dotted histogram 
represents the resuh for < < (1 - with M = 40 GeV. 



sented by the solid histogram for M = 40 GeV and by the dashed histogram for 
M = 60 GeV. By including the NLL order contributions, ambiguities due to the 
choice of the factorization scale become sufficiently small compared with the case at 
the LL order of QCD. An analytic calculation of the transverse momentum distribu- 
tion for one-gluon radiation cross section multiplied hy K = exp{as{Q'^) /2TrCp'7T'^) 
integrated over the range 40 GeV < Q < 60 GeV is also shown by the dotted curve. 
The generated cross section approximately reproduces the enhancement due to the 
next-to- next-to-leading order (NNLO) corrections {0{a1) terms) for the Drell-Yan 
cross section in the large qt region. 

Here, the transverse momenta of the initial-state partons with spacelike virtual- 
ity are generated using the effective two-body vertices, where the three-body decay 
functions are included as the boundary conditions, namely, Iq<P< /ns^'(~^^)i foi' 
the virtuality P of the outgoing partons in Eq. (2-13), with 

' Vqgg{zi,Z2,zy 
m{zuZ2,z)_ 

Here, the function L^^'^ is the coefficient of the mass singular term given by 



fm\^^'^2,z) = exp 



(3-2) 



Z1+Z2 \Zl + Z2 



(3-3) 



with zi + Z2 + z = 1, for the process q{ko) — )• g{h + h) + q{k) g{li) + g{l2) + q{k), 
where Pgg\z) is the LL order splitting function for the g ^ gg process, and the 
function Vqgg{zi, Z2-, z) is given in Ref. ITl]). Here, Zj (i = 1, 2) denotes the momentum 
fraction of parton momentum /j for the momentum /cq. 
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Fig. 3. Transverse momentum distributions for the virtual photon integrated over 40 GeV < Q < 
60 GeV with = 1800 GeV, where the NLL order contributions are taken into account. The 
sohd and dashed histograms represent the results for M = 40 GeV and M = 60 GeV obtained 
with the JC' scheme, respectively. The dotted curve represents the one-gluon radiation cross 
section multiplied hy K — exp{as{Q^)/2TvCFT^^) integrated over 40 GeV < Q < 60 GeV. 

The asymptotic form for the soft gluon radiation with the above choice gives 
f^s^ ~ zi for zi <C Z2,z^ which corresponds to the angular ordering condition, 
such as 6l-^^l^ <^o'2'^''^''^ small angle approximation, which is the case for 

§4. Summary and comments 

In this work, we have studied the Monte Carlo algorithm in hadron-hadron 
collisions, in which the initial-state radiation is generated using a parton shower 
model including the next-to-leading logarithmic (NLL) order contributions. 

We calculated hard scattering cross sections for the process qq — t- j*g with the 
jet calculus scheme, in which a kinematical constraint for the gluon radiation is taken 
into account. This method is called the JC' scheme. In this method, this kinematical 
constraint guarantees a proper phase space boundary for the subtraction terms. 

At the leading-logarithmic (LL) order approximation, we have shown that the 
transverse momentum distributions of the virtual photon depend on the factorization 
scale even the next-to-leading order (NLO) contribution in the hard scattering cross 
section is added in the large transverse region. We also pointed out that the trans- 

*■* Apart from the accuracy of the scaling violation generated by parton showers, the choice 
of the upper limit as P < [—k^)zi for zi <C Z2,z, which is a frame-independent form compared 
with the angular ordering condition, may be useful even at the LL order approximation in parton 
evolutions. In this choice, the ambiguity due to the kinematical boundary for the branching vertices 
is approximately eliminated. 



PS(NLL)+NLO M=40GEV 




± 
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verse momentum distributions depend on the kinematical boundary of the branching 
processes in the initial-state parton radiations, which cannot be determined at the 
LL order approximation. 

The kinematical conditions of the parton radiations and the remnant of the 
collinear region, in which the mass singularity is subtracted, are parts of the NLL 
order contributions. Therefore, these terms cannot be neglected at the NLL order 
accuracy. In the algorithm presented in this paper, these contributions are included 
in Monte Carlo calculations. 

Furthermore, the jet calculus scheme provides a matching between the initial- 
state radiations and the hard scattering cross sections. By including the NLL order 
contributions, we found that ambiguities due to the choice of the factorization scale 
become small compared with the case at the LL order of QCD. Furthermore, the 
kinematical boundary for the branching vertices, which is not determined at the LL 
order of QCD, is given by the three-body decay functions {0{a^) contributions for 
the branching vertex). 

For the large qt region, the cross section is enhanced, which is consistent with 
the K-factor ~ exp(as/27rCir7r^) for the Drell-Yan processes. The generated cross 
section approximately reproduces the results including the 0{a'l) contributions for 
the hard scattering cross sections.'^ 

The method presented in this paper may be useful for the construction of more 
accurate Monte Carlo algorithms, in which the parton radiation at the NLL order as 
well as the hard scattering cross sections at the NLO of QCD are taken into account. 

This method can be extended to the flavor singlet sector implementing the parton 
shower algorithm based on momentum distributions studied in Ref. [T0|) . 

It may be important to evaluate various processes by using Monte Carlo methods 
including the NLL order contributions with the NLO terms, such as VF+jet and 
tt+iet processes, which are backgrounds for new physics research as well as Higgs 
boson productions in hadron-hadron collisions. 
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Appendix A 

In Monte Carlo calculation of hard scattering contributions, we generate the 
photon virtuality and the transverse momentum qt with appropriate weight 
according to the hard scattering cross section. Although the hard scattering cross 
sections may be evaluated by using more accurate numerical methods, we implement 
a rather simple method for testing the algorithm. 

In order to evaluate the hard scattering cross section, we integrate Eq. (2-14) 
over the range bs < < as, with fixed s. Here, s is generated by the parton 
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showers. The energy scale of the running coupling constant in this region can be 
chosen as = fs. Therefore, the integrated cross section is given by 



;,[JC'](/) 
qq{NS) ' 



a,b) 



o. da 
df- 



[JC'](7) 
qq{NS) 



df 



dr—. Mrs, 0)K^_(^^)^(r). 



27r 



For the soft gluon region 1 — 2fM < f < 1, we consider two separate regions, 
one for 1 — tm < f < 1 (I = 51) and the other for 1 — 2fM < f < 1 — tm {I = S2). 
For I = SI, we have 



log(l - r) 
l-f 



(1-f) 



The factorized cross section has no M dependence, because the subtraction term is 
integrated over the range < — < s(l — f) due to the delta function in Eq. (2-10). 
Here, 



qqiNS) V / — -"-gg 

because the subtraction term is given by 



p[JC']{S) 



where K, 



[JC'KS) 
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t) and Fc 



[JC'](5) 



qq 



(e, f ) are given in Ref. [8]) for the singlet sector with 
infrared regularization so that the total momentum of partons in the initial-state is 
conserved. 

For I = S2, we have 



K[i?'M'^\T) = 2Cf 



qq{NS) 



LI -f 
l-f-TM 



2 l-f-TM 

log- 



[I + f) log 



1 - f 



since < s{l — f) in this region. Here, the subtraction term in Eq. (2-17) is 



where Fqq'^ ^^'^'*(e,f) is given in Ref. 



In order to evaluate the soft gluon contribution, we integrate over the range 
(1 — r]s) < f < 1, with fixed s. Here, r]s is a cutoff parameter satisfying tJs < 2fM 
with Tm = M'^/s. 

In this model, the soft gluon contributions, (/ = S), are resummed as 



[JC'](5), 
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with 



j,[JC'](S).. N 



exp 



.[JC'](S). 
"qqiNS) ^^''^ 



r]s)/ao{s,0) 



"qqiNS) 

The calculated results for the range 1 — "^s ^ f < 1 are given by 

JJC'KSI), 



.[JC'KS). 
%(iVS) v^' ^ 



a 



gq{NS) 



for < r]s < T]\i and 



[JC'](5) 
qq{NS) 



0" 



[JC'](51) 
gg-(iV5) 



(1, 1 - Tm) + ^, 



[JC'](52) 
qg(iV5) 



(1 - TM, 1 - 



for TM <r]s< '^TM- 

For To > 1 — 2fM(To = Qo/'^)' = s and = {Pq+Pq)T with the weight 

given by ao{s,0)R^^P^^^g^ [s, 1 — tq). Here, Qg is the minimum virtuality of photons 
to be generated. 

For Tq < 1 — 2fM, we compare a uniformly generated random number R{0 < 
R < 1) with the ratio R^^ defined by 



R(7 



[JC] ~ X 

%-(7^S)(l'^0) 



with 



+2(7, 



[JC'](C) 



qq 



(1 - 2fM,T0; 



Since (1— 2fM, 'Tq) ^ '''^^(Ars) (-'-"^'^A^'' ''^o); we can safely neglect the remnant 

of the hard collinear contribution cr^-^^^^ (1 — 2fM,'ro) in the factorized cross section, 

even at the NLL order accuracy.® On the other hand, in both the MS and MS' 
schemes, the remnant of the collinear region becomes negative. Therefore, strong 
cancellation may occur between the hard scattering contribution and the initial-state 
radiations .'^ 

If Ra > R, we set = s and qx = {pq + j%)r with the weight given by 

[JC] .X 

%-(^)(l'^o). 



If Ra < R,f is generated by solving the equation R = ^l^^jlg-ji'^, '^^^q^NS)^^'^^^ 
for given R with the weight c''qq(7V5)(-'-' '^o)- Here, the contribution of the hard gluon 
radiation (/ = H) integrated over the range ^ —t < s{l — f) — is given by 

^s(l-f)-M2 



JJC] 



d{-i)k![->{T. 



-i) 



2Cf 



(r) log '-^-^ - (1 - r - 2tm) 



Cf 



Tm 
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with 



k'(^{f^-i) = CF 



+ 



-i ' S{l-f)+ij Cf'"" 

The transverse momentum of the photon is constructed as gf. = {—t){s + t — 
Q'^)/s, where —i is generated according to the distribution of Kg^g^{f,—i) within 
M2 < -i< (1 -f)s -M2. 

Appendix B 



In this appendix, the exphcit expressions, which are implemented in the numer- 
ical calculations in the text with the JC' scheme, are presented. 

In order to simplify our analysis, we evaluate the cross section for the soft gluon 
radiation with the coupling constant a^lQ^) — The running coupling constant 

for the accuracy of the NLL order is normalized as ag (M|) = 0.114 at the Z° boson 
mass.'i^' 

The integrated cross section over the range {1 — r]s)s < Q'^ < s is given by 



JJC'](5) 
^qq{NS) 



(1, 1 - Vs 



2tt 



&o{s,0)I 



[JC'](5) 



qg{NS) 



(1,1 -r]s 



with 



Iqq{NS)(^^^-^s) 



^ ^^^[JC'](5) 

l-Vs 



The calculated results are given by 



JJC'](S) 



(1, 1-Vs)= I 



[Jc'Ksi) 

qq(NS) 



(1,1 



Cf 



45P_(r/„0) -21og2r/, + 31og 



- + -tt' 
2 3 



for < rjs < tm and 



r[JC'](5) 
^qq(NS) ' 



1, ^-Vs)= I 



[JC'](51) 
qq(NS) ' 



1, 1 - tm) + / 



[JC'](52) 
qq{NS) ' 



1 - TM, l-ris 



with 



i^iZE^\i-fM,i-%) = 2CF 



log 



1 - TM 



log 



2SP- 



lis - TM 
1 - TM 



,0 +25P4 



Vs - TM 
TM 



,0 -5P_fe,fM) 



(1 - tmY 



+ 1 - -TA/ + 2(7?s - Tm) + Vs lof 



TM 



for TM <r]s< '^TM- 

Here, we define 



SP±{a,b) 
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The contribution of the hard gkion radiation (/ = H) integrated over the range 
< -i < s{l - f ) - and tq < r < 1 - Itm is given by 

cr^g- '(1 - 2tm, To) ~ — ao{s, 0)/^^- '(1 - 2tm,tq) 



with 



4f (l-2fM,ro) = 2C7^ 



TM 



1 - TM ' '^-TM 



+ log^^flogl^-l + 2fMV(l-fM--fo)log ^-t^-^-"° +2(l-2fM-fo)' 



To 



TM 



with an averaged couphng constant as(s,a, 6) = (as(as) + as(6s))/2, where a and 6 
are the upper and lower hmits of f integration, respectively. 



2 
3 

4' 
5 
6 

r 
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